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1.  Introduction 


Micro-power  generation  technologies  have  recently  received  a  great  deal  of  attention  from  both 
government  laboratories  and  academia.  The  majority  of  the  work  under  investigation  today  has 
focused  on  the  miniaturization  of  macro-scale,  fuel  dependant  power  technologies  such  as  micro¬ 
turbines  (1 )  and  micro-generators,  (2).  Micro-power  technologies  that  use  their  environment  for 
their  sources  of  input  energy  (3-6),  as  opposed  to  requiring  on  board  fuel,  should  provide  the 
benefits  of  longer  operational  lifetimes  and  significantly  lower  system  weight  and  volume.  One 
detriment  to  these  other  approaches  is  the  energy  conversion  loss  associated  with  transformations 
from  multiple  energy  domains.  Microsystems  that  combust  fuels  (7),  for  example,  transform 
energy  from  chemical  to  thermal  to  mechanical  to  electrical  domains;  the  second  law  of 
thermodynamics  dictates  significant  energy  transformation  losses  that  necessitate  greater 
quantities  of  fuel,  thus  increasing  the  system  size  and  volume  and  negatively  impacting 
operational  lifetimes.  Piezoelectric  vibrational  energy  reclamation  requires  only  one  energy 
transformation. 

This  work  demonstrates  the  feasibility  of  vibrational  energy  harvesting  via  thin  film  PZT 
MEMS.  Using  the  direct  piezoelectric  effect,  periodic  deflection  of  a  piezoelectric  thin  film 
produces  a  generated  charge  that  can  be,  for  each  deflection  cycle,  stored  for  later  external 
system  power  requirements.  The  input  energy  is  reclaimed  from  the  waste  vibrational  energy  of 
another  mechanical  system  vibrating  at  a  given  frequency  and  amplitude.  PZT  is  the 
piezoelectric  thin  film  of  choice  for  this  application,  given  its  superior  electro-mechanical 
coupling  coefficient  (8). 


Figure  1.  Prototype  piezoelectric  vibrational  energy 
reclamation  device. 

We  present  analytical  device  modeling,  MEMS  fabrication,  and  the  results  from  laser  Doppler 
vibrometry  and  vibration  testing.  Piezoelectric  charge  generation  derived  from  piezoelectric 
constitutive  equations  (9),  a  discretized  Bernoulli-Euler  mechanical  deformation  model,  and 
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simple  power  density  derivation  provide  the  basis  for  the  modeling  effort.  Laser  Doppler 
vibrometry  (10)  allows  for  the  characterization  of  the  device  frequency  response  and  resonant 
mode  shape  investigation,  while  the  use  of  an  electro-mechanical  vibration  table  enables  direct 
measurement  of  the  piezoelectrically  generated  current.  The  basic  geometry,  coordinate  system, 
and  composition  of  the  cantilever  device  are  shown  in  figure  2. 


Figure  2.  Conceptual  cross-section  of  cantilever 
vibrational  energy  harvester,  displaying 
coordinate  system  employed. 


2.  Vibrational  Energy  Reclamation  Design 


2.1  Piezoelectric  Charge  Generation 

The  charge  generated  by  a  piezoelectric  cantilever  beam  is  obtained  by  integrating  the  electric 
displacement  perpendicular  to  the  electrodes  D3  over  the  area  of  the  beam.  The  following  is  an 
approximation  of  D3: 

D3  =  £?31^Wxx  (1) 

t?piezo  =  |D3h(x)dx  (2) 

<7Piezo  =  e3i/?jwxxb(x)dx  (3) 

The  approximation  arises  from  the  presence  of  the  h  term,  the  distance  from  the  neutral  axis 
(N.A.)  to  the  midplane  of  the  piezoelectric  layer.  These  expressions  also  assume  uniform 
deformation  across  the  width  of  the  beam.  By  employing  Bemoulli-Euler  beam  theory,  the  strain 
distribution  is  assumed  linear  with  respect  to  the  device  thickness;  implying  the  average  strain  is 
referenced  by  the  h  term.  For  a  constant  width  electrode,  equation  3  then  becomes: 

<7  piezo(t)  =  e.j\  bhQn(l)  (4) 
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2.2  Bernoulli-Euler  Mechanical  Deformation  Model 


The  devices  are  electrically  coupled  piezoelectric  cantilever  arrays.  Modeling  the  basic 
cantilever  structure  as  a  Bernoulli-Euler  beam  pennits  the  beam  to  be  discretized  into  an 
arbitrary  number  of  elements  (FEM)  with  the  deflections  at  each  node  described  by,  (11): 


(5) 


where  the  vector  on  the  left  of  equation  5  contains  the  vertical  and  angular  displacements  of  each 
node  (the  last  tenn  is  the  angular  deflection  required  to  solve  the  charge  equation  4).  The  vector 
on  the  right  contains  the  amplitudes  of  the  forces  and  moments  at  each  node.  Both  vectors  are 
(2n  x  1)  vectors  where  n  is  the  number  of  nodes  the  beam  has  been  discretized  into.  [M]  and  [K], 
the  global  mass  and  stiffness  matrices,  are  assembled  from  their  respective  elemental  matrices. 
For  the  one-dimensional  case,  the  elemental  matrices  are  assembled  along  the  diagonal  of  an 
empty  (2n  x  2n)  matrix,  with  adjoining  nodal  values  summing.  Equations  6  and  7  are  the  one¬ 
dimensional  elemental  mass  and  stiffness  matrices;  assuming  a  cubic  shape  function  for  the 
elements,  (11). 
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2.2.1  Conversion  of  Base  Motion  to  Uniformly  Distributed  Inertial  Harmonic  Load 

A  cantilever  beam  subjected  to  hannonic  motion  of  its  base  is  acted  upon  by  a  uniform 
distributed  inertial  harmonic  excitation  force  equal  to  equation  8,  (12): 

F  =  7cQ2mbeamSin(cot)  (8) 

In  order  to  represent  this  uniform  distributed  inertial  load  in  the  discretized  Bernoulli-Euler 
equation  of  motion,  the  force  must  be  adequately  discretized  along  the  length  of  the  beam.  In 
motion  equation  5,  the  elemental  force  tenns  in  the  force  and  moment  vector  become  the 
amplitude  of  F  divided  by  the  number  of  elements. 
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2.2.2  Determination  of  Composite  Beam  Modulus  of  Elasticity,  Moment  of  Inertia,  Unit 
Mass,  Neutral  Axis  Location,  and  Resonance 

The  area  moment  of  inertia  of  a  cross-section  of  the  beam  with  respect  to  the  y-axis  of  the 
composite  beam  is  the  summation  of  the  sums  of  the  moments  of  inertia  of  each  layer  and  the 
product  of  each  layer's  cross-sectional  area  and  the  square  of  the  distance  from  the  neutral  axis  to 
the  centroidal  axis  of  each  layer,  as  given  by  the  parallel  axis  theorem  (13). 

2 

-^composite  =  L(Il+AldO  (9) 

The  composite  modulus  of  elasticity  is  the  summation  of  the  products  of  volume  fractions  of 
each  layer  and  the  modulus  of  the  corresponding  layer.  For  a  constant  width  beam,  this  reduces 
to  (13): 


^composite  (l layerT'layer  9 beam) 


(10) 


The  location  of  the  neutral  axis  is  found  by  dividing  the  summation  of  the  products  of  the  layer 
cross  sectional  areas  and  the  distances  from  the  layer  centroids  to  some  arbitrary  reference  axis 
by  the  cross-sectional  area  of  the  entire  beam,  (13). 


Vbar  (E  (^4iVi)V^ltotal 


(11) 


The  mass  per  unit  length  is  found  by  multiplying  the  beam  width  by  the  summation  of  the 
products  of  the  layer  densities  and  the  layer  thickness. 


m beam  b\  E  ( /t a\  e;t la\ or  )]  (12) 

The  natural  frequencies  of  the  beams  are  the  eigenvalues  of  the  following  matrix  expression: 

[M]-l[K]{q0}  =  an2{q0}  (13) 


Figure  3.  Micrograph  of  released  cantilever  test  structures. 
Residual  film  stress  deformation  is  apparent. 
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2.3  Power  Density 

The  average  power  developed  by  a  piezoelectric  cantilever  structure  is  equal  to  half  of  the 
maximum  voltage  across  the  PZT  multiplied  by  the  maximum  current  produced  by  the  device. 
The  voltage  amplitude  equation  16,  across  the  PZT,  is  equal  to  the  generated  charge  divided  by 
the  capacitance  equation  14  of  the  PZT  layer.  Assuming  that  the  charge  is  completely  extracted 
once  during  the  deflection  cycle,  the  current  amplitude  equation  15  is  the  product  of  the 
operating  frequency  in  hertz  and  the  charge  generated  per  deflection  cycle. 


Cp  (£33  £0  l  b)/tc 

(14) 

ipiezo(t)  =  q  cocos(cot) 

(15) 

hpiezo(t)  =  ((<?31  h  0n  4)/(£ 33  £0  /)) sin(cot) 

(16) 

Power  =  (q  ipiezo)/2Cp 

(17) 

Power/ Area  =  (e-n  h  0n)2  co  tc! 2e33  £0  /2 

(18) 

2.4  Mitigation  of  Residual  Stress  Deformation 

In  order  to  investigate  the  behavior  of  the  longer  MEMS  devices,  on  the  order  of  a  few 
millimeters,  it  was  necessary  to  compensate  for  the  residual  thin  film  stress  deformation  of  these 
highly  compliant  structures.  If  the  mean  force  associated  with  the  stress  gradient,  whose  location 
can  be  detennined  from  the  mean  value  theorem  of  integrals,  acts  at  a  distance  displaced  from 
the  structure’s  N.A.,  residual  stress  deformation  will  occur.  This  stress  gradient,  through  the 
thickness  of  the  device,  leads  to  the  distributed  moment  that  serves  to  deform  the  structure  out  of 
plane.  This  deformation  can  be  so  extreme  as  to  prevent  either  measurement  of  a  device  or  even 
device  operation.  By  designing  the  structure’s  film  thicknesses  and  engineering  the  stress 
gradients,  a  near  planar  structure  can  be  achieved. 


3.  Fabrication 


We  fabricated  PZT-based  cantilevers  and  cantilever  arrays  using  a  simple,  four  mask  fabrication 
process.  A  schematic  cross-section  of  a  completed  device  is  illustrated  in  figure  2.  Starting  with 
a  bare  Si  wafer,  a  micron  and  a  half  of  PECVD  SiC>2  was  deposited  at  250  °C.  The  film  was  then 
annealed  at  700  °C  in  a  nitrogen  atmosphere.  An  adhesion/diffusion  barrier  layer  of  titanium 
was  sputtered  on  the  PECVD  oxide,  followed  by  a  sputtered  platinum  bottom  electrode  layer. 
Sol-gel  PZ0.52T0.48  (PZT)  was  then  repeatedly  spun,  pyrolized,  and  crystallized  until  the 
appropriate  PZT  thickness  was  achieved.  A  liftoff  of  sputtered  platinum  defined  the  top 
electrodes  and  bond  pads. 
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Table  1 .  Thin  film  thickness  for  the  devices 
tested;  measured  via  Variable  Angle 
Ellipsometry. 


Material 

Thickness  (nmj 

Silicon  Dioxide 

1688 

Titanium 

20 

Bottom  Platinum 

120 

PZT 

2048 

Top  Platinum 

100 

After  photo  patterning,  we  defined  the  release  trench  and  bottom  electrode  region  by  ion  milling 
of  the  PZT.  The  remaining  photoresist  was  ashed  in  an  oxygen  plasma.  The  now  exposed 
platinum  in  the  release  trench  and  an  electrical  isolation  trench  were  ion  milled  through  the 
platinum  and  underlying  Si02,  exposing  the  underlying  silicon  substrate.  After  ashing,  a  fourth 
photolithography  step  defining  the  XeF2  release  windows  was  perfonned.  Given  the  1:1  etch 
aspect  ratio  of  the  XeF2,  the  release  windows  were  positioned  a  distance  Vi  the  beam  width  from 
the  anchors.  This  insured  that  the  actual  beam  lengths  corresponded  to  the  designed  lengths 
when  the  beams  released.  Device  die  were  cleaved  and  released  individually,  undercutting  the 
beams.  After  releasing  in  XeF2,  the  photoresist  was  again  ashed. 


4.  Experimental  Results 


Multiple  devices  and  cantilever  test  structures  with  resonant  frequencies  from  570  Hz  to  240  kHz 
were  fabricated  using  the  process  flow  described  in  section  3.  The  fabricated  devices  and  test 
structures  were  packaged  in  forty  pin  dual  inline  package  (DIPs),  compatible  with  both  testing 
apparatuses  employed  in  device  characterization,  laser  Doppler  vibrometry  and  electro¬ 
mechanical  vibration  excitation.  Figures  4(a)  and  4(b)  display  the  experimental  setups  used  for 
this  work. 


Figure  4.  Block  diagrams  of  the  experimental  setups  used  in  characterizing  piezoelectric  vibrational 
energy  reclamation  devices;  (a)  laser  Doppler  vibrometer  and  (b)  the  vibration  table. 
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4.1  LDV  Frequency  Response/Mode  Shape  Evaluation 

The  packaged  devices  were  first  characterized  in  terms  of  their  frequency  response.  The  devices 
were  poled  at  their  coercive  field  for  approximately  one  hour.  They  were  excited 
piezoelectrically  at  1  Vp.p  and  measured  with  a  Polytec™  laser  Doppler  vibrometer  at  atmospheric 
pressure.  A  detailed  discussion  of  this  measurement  technique  can  be  found  in  other  references 
(10). 


Figure  5  and  table  2  demonstrate  good  agreement  between  the  eigen  solutions  of  equation  13  and 
the  LDV  experimentally  verified  fundamental  resonant  frequencies.  The  underestimating  of  the 
natural  frequencies  may  be  due  to  a  number  of  issues,  including  errors  in  the  Young’s  modulus 
terms  used  in  equation  13,  the  presence  of  residual  film  stresses  (14),  air  damping,  and  deviation 
of  actual  beam  length  from  designed/assumed  lengths  associated  with  the  release  process.  Even 
with  these  uncertainties,  the  values  are  still  within  a  few  percent  of  experimental  values. 
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Figure  5.  Experimental  and  predicted  fundamental 
resonant  frequency  for  varied  cantilever 
beam  lengths. 


Table  2.  Error  between  predicted  and  experimental  fundamental 
natural  frequencies. 


Beam  Length  (Microns) 

%  Error  in  Natural  Freq. 

400 

0.7 

500 

1.2 

600 

3.6 

700 

1.6 

800 

1.2 

900 

2.4 

1000 

1.25 

The  frequency  response,  up  to  100  kHz,  of  a  2  mm  long  cantilever  and  an  1100  pm  long 
cantilever  are  shown  in  figures  6a  and  6b.  The  shift  in  predicted  response  is  more  apparent  at 
higher  frequencies,  as  the  error  (table  2)  through  higher  order  modes  remains  roughly  constant. 
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In  addition  to  frequency  response  behavior,  laser  Doppler  vibrometry  allows  measurement  of  the 
resonant  mode  shapes.  Figure  7  displays  the  fundamental  harmonic,  7th  order  mode  shapes 
(region  near  the  anchor  only  displayed  due  to  scan  area  limitations)  of  a  2  mm  long  device, 
respectively.  The  observed  mode  shapes  agree  well  with  expectations.  However,  torsional 
modes  of  resonance  (figure  7)  were  also  observed.  Their  occurrence  does  appear  to  alter  the 
measured  frequency  response,  as  the  current  model  does  not  account  for  torsional  modes.  For 
the  devices  tested,  torsional  resonant  modes  are  present  in  the  range  tested  with  the  vibration 
table,  as  they  were  identified  with  the  Laser  Doppler  Vibrometer.  A  more  detailed  investigation 
of  these  modes  and  their  impact  on  the  frequency  response  shall  be  pursued  in  future  work. 


Figure  6.  Experimental  and  predicted  frequency  responses  of  a  2  mm  long  (a)  and  an  1100  pm  long  device  (b); 
the  top  plot  in  each  is  experimental  data  obtained  from  LDV  while  the  lower  plot  is  the  predicted 
response  w/arbitrary  magnitudes. 


Figure  7.  Fundamental,  7th  order,  and  torsional  resonant  mode  shapes  of  a  2  mm  long  cantilever  (only  a  portion  near 
the  anchor  is  actually  visible). 

4.2  Vibration  Table  Analysis 

A  B&K  type  4808  vibration  table,  driven  by  a  type  2712  power  amplifier,  was  used  to 
mechanically  excite  the  piezoelectric  cantilever  test  structure  packages.  The  package  was 
mounted  in  a  specially  designed  fixture,  containing  a  forty-pin  zero  insertion  force  (ZIF)  socket, 
which  allowed  multiple  packages  to  be  tested  with  minimal  work.  An  accelerometer  was 
mounted  to  the  fixture  to  measure  the  vibrational  amplitudes.  The  table’s  fundamental  resonance 
was  10.1  kHz,  well  above  the  range  investigated.  Data  was  taken  between  100  Hz  and  1  kHz, 
and  at  any  resonance  peaks  present  in  the  range.  A  lock-in  amplifier  was  used  to  drive  the  table 
and  to  measure  the  generated  piezoelectric  current  and  voltage  from  the  devices.  As  the 
amplitude  of  vibration  was  not  easily  decoupled  from  the  frequency,  the  data  displayed  in 
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figure  8,  represents  variable  amplitudes  of  vibration.  Figures  8a,  8b,  and  8c  display 
piezoelectrically  generated  current  from  vibration  table  testing  for  an  1100  pm  long,  1500  pm 
long,  and  2  mm  cantilevers,  respectively.  There  is  reasonable  agreement  between  the  predicted 
and  experimental  currents.  The  deviation  is  thought  to  be  due  in  part  to  the  slight  errors  in 
resonant  frequencies  and  the  presence  of  torsional  modes  contributing  to  the  frequency  response. 
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Figure  8.  Piezoelectrically  generated  current  for  an  1 100  pm  long  device  (a),  a  1500  pm  long  device  (b),  and  a 

2  mm  long  device  (c)  from  vibrational  excitation;  model  vs.  experimental.  Note:  Displacement  is  coupled 
to  frequency  in  the  data. 


The  power  generated  at  570  Hz  at  a  peak-to-peak  displacement  of  3.4  pm  by  the  2  mm  long 
device  was  0.16  nW.  This  translates  into  0.08  pW/ctn  .  The  Laser  Doppler  Vibrometer  and 
mechanical  excitation  testing  of  the  prototype  devices  have  demonstrated  the  validity  of  the 
current  model.  The  model  predicts  power  generation  of  1.54  mW/cm"  at  700  Hz/50  pm 
amplitude  of  vibration  for  an  optimized  device. 


5.  Conclusions 


Surface  micromachined  piezoelectric  PZT  vibrational  energy  reclamation  devices  have  been 
designed,  fabricated,  and  tested  to  demonstrate  the  feasibility  of  vibrational  energy  reclamation 
via  thin  fdm  PZT  MEMS.  Experimental  data  has  validated  the  developed  electro-mechanical 
model.  Predictions  of  this  model  imply  a  power  density  of  1.54  mW/cnT  at  a  700  Hz,  50  pm 
amplitude  vibration.  Future  work  will  focus  on  model  optimization,  development  of  an 
appropriate  charge  extraction  circuit,  and  system  demonstration. 
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Electric  displacement  perpendicular  to  the  electrodes 
Piezoelectrically  generated  charge 
Neutral  Axis 

Piezoelectric  stress  constant 

Distance  from  the  N.A.  to  the  PZT  midplane 

Beam  width 

Beam  curvature,  2nd  derivative  of  deflection 
Beam  Axis  dimension 

Angular  deflection  at  the  free  end  of  the  cantilever  beam 

Global  mass  matrix 

Global  stiffness  matrix 

Operating  frequency  in  hertz 

Amplitude  of  vibration 

Displacement  Vector 

Mass  of  composite  beam 

Total  inertial  force 

Location  of  N.A.  relative  to  arbitrary  reference  axis 
Young’s  Modulus 
Moment  of  Inertia 
Mass  density 

Thickness  of  piezoelectric  layer 

Piezoelectric  layer  dielectric  constant 

Pennittivity  of  free  space 
(8.854  x  1CT12  F/m) 

Piezoelectrically  generated  current 
Length  of  cantilever 
Piezoelectric  layer  capacitance 
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